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1 Introduction 



One of the most important motivations for studies of D° — D° mixing is the possibility of 
observing a signal from new physics which can be separated from the one generated by the 
Standard Model (SM) interactions. The D° - D° mixing proceeds extremely slowly, which 
in the Standard Model is usually attributed to the absence of superheavy quarks destroying 
GIM cancelations Q . The low energy effect of new physics particles can be naturally written 
■ in terms of a series of local operators of increasing dimension generating AC = 2 transitions. 

. These operators, as well as the one loop Standard Model effects, generate contributions to 

the effective operators that change D° state into D° state leading to the mass eigenstates 

^ \D 1 )=p\D°)± q \D°), (1) 

Ph! _ 

O ■ where the complex parameters p and q are obtained from diagonalizing the D — D° mass 

matrix. The mass and width splittings between these eigenstates are parameterized by 

x = {rn 2 - mi )/T, y = {Y 2 -T x )/{2Y), (2) 

•i— I , 

^ ■ where and Ti^ are the masses and widths of D± t 2 and the mean width and mass are 

r = (Ti + T2)/2 and m = (mi + m,2)/2. Since y is constructed from the decays of D into 
physical states, it should be dominated by the Standard Model contributions, unless new 
physics significantly modifies AC = 1 interactions. 

Presently, experimental information about the D° — D° mixing parameters x and y 
comes from the time-dependent analyses that can roughly be divided into two categories. 
First, more traditional studies look at the time dependence of D — > f decays, where / is 
the final state that can be used to tag the flavor of the decayed meson. The most popular is 
the non-leptonic doubly Cabibbo suppressed decay (DCSD) D° — ► K + ir~ . Time-dependent 
studies allow one to separate the DCSD from the mixing contribution D° — > D° — > K + ir~ , 



r[D (t)^K+ir-} = e- rt \A K - w+ \ 2 

R2 _ J 

(3) 



i? 2 

R + VRR m (y' cos0 - x' sin <t>)Tt + ~^{y 2 + x 2 ){Tt) 2 
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where R is the ratio of DCS and Cabibbo favored (CF) decay rates and q/p = R m e 1 ^. Since 
x and y are small, the best constraint comes from the linear terms in t that are also linear in 
x and y. A direct extraction of x and y from Eq. (|3|) is not possible due to unknown relative 
strong phase 5 of DCS and CF amplitudes Q , as x' = x cos 5 + y sin 5, y' = y cos 5 — x sin 5. 
This phase can be measured independently [[|. The corresponding formula can also be 
written for D° decay with x' — ► —x' and R m — > i?" 1 Q. 

Second, Z?° mixing can be measured by comparing the lifetimes extracted from the 
analysis of D decays into the CP-even and CP-odd final states. This study is also sensitive 
to a linear function of y via 



t(D -> K~ir + ) 
t(D -> K+K~) 



1 



y cos i 



x sin ( 



R~ 
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(4) 



Time-integrated studies of the semileptonic transitions are sensitive to the quadratic form 
x 2 + y 2 and at the moment are not competitive with the analyses discussed above. The 
construction of a new tau-charm factory at Cornell will introduce other time-independent 
methods that are sensitive to a linear function of y ||. 

The current experimental upper bounds on x and y are on the order of a few times 



1(T 



and are expected to improve significantly in the coming years. To regard a future 



discovery of nonzero x or y as a signal for new physics, we would need high confidence that 
the Standard Model predictions lie well below the present limits. As was recently shown 0, 
in the Standard Model x and y are generated only at second order in SU(3) breaking, 



x , y ~ sin 2 C x [SU{3) breaking]' 



(5) 



where 9c is the Cabibbo angle. Therefore, predicting the Standard Model values of x and 
y depends crucially on estimating the size of SU(3) breaking. Although y is expected to 
be determined by Standard Model processes, its value nevertheless affects significantly the 
sensitivity to new physics of experimental analyses of D mixing 0]. 



2 Theoretical Expectations 

Theoretical predictions of x and y within and beyond the Standard Model span several 
orders of magnitude (see Ref. 0). Roughly, there are two approaches, neither of which 
give very reliable results because m c is in some sense intermediate between heavy and 
light. The "inclusive" approach is based on the operator product expansion (OPE). In the 
m c 3> A limit, where A is a scale characteristic of the strong interactions, AM and Ar 
can be expanded in terms of matrix elements of local operators || . Such calculations yield 
x,y < 1CT 3 . The use of the OPE relies on local quark-hadron duality, and on A/m c being 
small enough to allow a truncation of the series after the first few terms. The charm mass 
may not be large enough for these to be good approximations, especially for nonleptonic D 
decays. An observation of y of order 10~ 2 could be ascribed to a breakdown of the OPE or 
of duality, but such a large value of y is certainly not a generic prediction of OPE analyses. 
The "exclusive" approach sums over intermediate hadronic states, which may be modeled 



2 



or fit to experimental data [|]. Since there are cancellations between states within a given 
SU(3) multiplet, one needs to know the contribution of each state with high precision. 
However, the D is not light enough that its decays are dominated by a few final states. In 
the absence of sufficiently precise data on many decay rates and on strong phases, one is 
forced to use some assumptions. While most studies find x,y < 1CP 3 , Refs. |J obtain x 
and y at the 10 -2 level by arguing that SU (3) violation is of order unity, but the source of 
the large SU(3) breaking is not made explicit. 

In what follows we first prove that D° — D° mixing arises only at second order in SU(3) 
breaking effects. The proof is valid when SU(3) violation enters perturbatively. This would 
not be so, for example, if D transitions were dominated by a single narrow resonance close to 
threshold @, ffi|. Then we argue that reorganization of "exclusive" calculation by explicitly 
building SU(3) cancellations into the analysis naturally leads to values of y ~ 1% if only 
one source of SU(3) breaking (phase space) is taken into account. 

The quantities M\i and Y\i which determine x and y depend on matrix elements 
(D°\ TL w TLw \D°) , where H w denote the AC = — 1 part of the weak Hamiltonian. Let 
D be the field operator that creates a D° meson and annihilates a D°. Then the matrix 
element, whose SU(3) flavor group theory properties we will study, may be written as 

{0\DH w H w D\0). (6) 

Since the operator D is of the form cu, it transforms in the fundamental representation of 
SU(3), which we will represent with a lower index, D{. We use a convention in which the 
correspondence between matrix indices and quark flavors is (1,2,3) = (u,d,s). The only 
nonzero element of D{ is D\ = 1. The AC = — 1 part of the weak Hamiltonian has the 
flavor structure (jlic)(qjqk), so its matrix representation is written with a fundamental index 
and two antifundamentals, HV . This operator is a sum of irreps contained in the product 
3x3x3 = 15 + 6 + 3 + 3. In the limit in which the third generation is neglected, HV is 
traceless, so only the 15 and 6 representations appear. That is, the AC = — 1 part of TC W 
may be decomposed as \{Oj^ + Oq), where 

Oj£ = (sc)(ud) + (uc)(sd) + si(dc)(ud) + si(uc)(dd) 

— si(sc)(us) — s\(uc)(ss) — s1(dc)(us) — sf(uc)(ds) , 
06 = (sc)(ud) — (uc)(sd) + si(dc)(ud) — si(uc)(dd) 

— si(sc)(us) + si(uc)(ss) — sf(dc)(us) + s1(uc)(ds) , (7) 

and si = sin 9c- The matrix representations H(15)^ and i?(6)Y have nonzero elements 



H (I5)« : Hf = Hf = 1 , H? = Hf = s x , 

Hf = Hf = -s x , HP = Hf = -si , 



HP = -HP = - Sl , = -Hf = -si . 

We introduce SU(3) breaking through the quark mass operator A4, whose matrix repre- 
sentation is Mj = diag(m u , md,m s ) as being in the adjoint representation to induce SU(3) 
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violating effects. We set m u = nid = and let m s ^ be the only SU(3) violating param- 
eter. All nonzero matrix elements built out of D{, and Mj must be SU(3) singlets. 

We now prove that D° — D° mixing arises only at second order in SU(3) violation, by 
which we mean second order in m s . First, we note that the pair of D operators is symmetric, 
and so the product D{Dj transforms as a 6 under SU(3). Second, the pair of Ti w 's is also 
symmetric, and the product H^H ™ is in one of the reps which appears in the product 

[(15 + 6) x (l5 + 6)] 5 = (15 x l5) 5 + (15 x 6) + (6 x 6)s (9) 

= (60 + 24 + 15 + 15' + 6) + (42 + 24 + 15 + 6 + 3) + (15' + 6) . 

A direct computation shows that only three of these representations actually appear in the 
decomposition of Ti w 7i w . They are the 60, the 42, and the 15' (actually twice, but with 
the same nonzero elements both times). So we have product operators of the form (the 
subscript denotes the representation of SU(3)) 

dd = v 6 , n w n w = o m + 042 + Ow • (10) 

Since there is no 6 in the decomposition of 7i w 7i w , there is no SU(3) singlet which can be 
made with 2?g, and no SU(3) invariant matrix element of the form (JBJ) can be formed. This 
is the well known result that D°—D° mixing is prohibited by SU(3) symmetry. Now consider 
a single insertion of the SU{3) violating spurion M. The combination DqM transforms as 
6 x 8 = 24 + 15 + 6 + 3. There is still no invariant to be made with Ti. w 7i w , thus D° — D° 
mixing is not induced at first order in SU (3) breaking. With two insertions of M, it becomes 
possible to make an SU(3) invariant. The decomposition of VMM is 

6 x (8 x 8)5 = 6 x (27 + 8 + 1) = (60 + 42 + 24 + 15 + TE' + 6) + (24 + 15 + 6 + 3) + 6 . (11) 

There are three elements of the 6 x 27 part which can give invariants with 7i w 7i w . Each 
invariant yields a contribution to D° — D° mixing proportional to s\m1. Thus, D° — D° 
mixing arises only at second order in the SU(3) violating parameter m s . 

We now turn to the contributions to y from on-shell final states, which result from every 
common decay product of D° and D°. In the SU(3) limit, these contributions cancel when 
one sums over complete SU(3) multiplets in the final state. The cancellations depend on 
SU(3) symmetry both in the decay matrix elements and in the final state phase space. 
While there are SU(3) violating corrections to both of these, it is difficult to compute the 
SU(3) violation in the matrix elements in a model independent manner. Yet, with some mild 
assumptions about the momentum dependence of the matrix elements, the SU(3) violation 
in the phase space depends only on the final particle masses and can be computed. We 
estimate the contributions to y solely from SU(3) violation in the phase space. We find 
that this source of SU(3) violation can generate y of the order of a few percent. 

The mixing parameter y may be written in terms of the matrix elements for common 
final states for D° and D° decays, 

v = f E / P- s -]« <^°l n - n ™ 1^°) > ( 12 ) 

n J 
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where the sum is over distinct final states n and the integral is over the phase space for 
state n. Let us now perform the phase space integrals and restrict the sum to final states 
F which transform within a single SU(3) multiplet R. The result is a contribution to y of 
the form 

(D°\H w Lcp(Fr) £ \n)pn(n\\n w \D°), (13) 

*- n£F R > 

where p n is the phase space available to the state n, r/cp = ±1 §§. In the 57/(3) limit, all 
the p n are the same for n G Fr, and the quantity in braces above is an SU(3) singlet. Since 
the p n depend only on the known masses of the particles in the state n, incorporating the 
true values of p n in the sum is a calculable source of SU (3) breaking. 

This method does not lead directly to a calculable contribution to y, because the matrix 
elements {n\H w \D Q ) and (D°\H w \n) are not known. However, CP symmetry, which in the 
Standard Model and almost all scenarios of new physics is to an excellent approximation 
conserved in D decays, relates (D°\H w \n) to (D°\7i w \n). Since \n) and |n) are in a common 
SU(3) multiplet, they are determined by a single effective Hamiltonian. Hence the ratio 

_ En£F R (D°\H W \n)pn(n\H w \D ) = En e F R (D°\H W \n)pn(n\H w \D ) 
VF ' R ZneF R (D°\H w \n)p n (n\H w \D°) E^rpo - n ) ( ' 

is calculable, and represents the value which y would take if elements of Fr were the only 
channel open for D decay. To get a true contribution to y, one must scale yF,R to the 
total branching ratio to all the states in Fr. This is not trivial, since a given physical final 
state typically decomposes into a sum over more than one multiplet Fr. The numerator of 
yF,R is °f order sf while the denominator is of order 1, so with large SU(3) breaking in the 
phase space the natural size of yF,R is 5%. Indeed, there are other 577(3) violating effects, 
such as in matrix elements and final state interaction phases. Here we assume that there 
is no cancellation with other sources of SU(3) breaking, or between the various multiplets 
which occur in D decay, that would reduce our result for y by an order of magnitude. This 
is equivalent to assuming that the D meson is not heavy enough for duality to enforce such 
cancellations. Performing the computations of yF,R H> we see that effects at the level of 
a few percent are quite generic. Our results are summarized in Table |l|. Then, y can be 
formally constructed from the individual yF,R by weighting them by their D° branching 
ratios, 

1 



V 



y = ^ 2^ yp,R 



(15) 

F,R L neF fl 

However, the data on D decays are neither abundant nor precise enough to disentangle the 
decays to the various 577(3) multiplets, especially for the three- and four-body final states. 
Nor have we computed yF,R for all or even most of the available representations. Instead, 
we can only estimate individual contributions to y by assuming that the representations 
for which we know yFR to be typical for final states with a given multiplicity, and then 
to scale to the total branching ratio to those final states. The total branching ratios of 
D° to two-, three- and four-body final states can be extracted from Ref. [11]. Rounding 



to the nearest 5% to emphasize the uncertainties in these numbers, we conclude that the 



5 



Final state representation 


Vf,r/si 


VF,R (%) 


PP 


8 


-0.0038 


-0.018 




27 


-0.00071 


-0.0034 


PV 


8a 


0.032 


0.15 




85 


0.031 


0.15 




10 


0.020 


0.10 




TO 


0.016 


0.08 




27 


0.04 


0.19 


(VV)s-wave 


8 


-0.081 


-0.39 




27 


-0.061 


-0.30 


(VV)p-wave 


8 


-0.10 


-0.48 




27 


-0.14 


-0.70 


( ^ ^ ' c/- wave 


8 


0.51 


2.5 




27 


0.57 


2.8 


(3P)s-wave 


8 


-0.48 


-2.3 




27 


-0.11 


-0.54 


(3P)p-wave 


8 


-1.13 


-5.5 




27 


-0.07 


-0.36 


( ^ -P ) f o r m . f act or 


8 


-0.44 


-2.1 




27 


-0.13 


-0.64 


4P 


8 


3.3 


16 




27 


2.2 


11 




27' 


1.9 


9.2 



Table 1: Values of uf,r for some two-, three-, and four-body final states. 



branching fractions for PP, (VV)s-wavei (^^)<i-wave anc ^ approximately amount to 
5%, while the branching ratios for PV and 4P are of the order of 10% ||. 

We observe that there are terms in Eq. (|l^), like nonresonant 4P, which could make 
contributions to y at the level of a percent or larger. There, the rest masses of the final state 
particles take up most of the available energy, so phase space differences are very important. 
One can see that y on the order of a few percent is completely natural, and that anything 
an order of magnitude smaller would require significant cancellations which do not appear 
naturally in this framework. Cancellations would be expected only if they were enforced 
by the OPE, or if the charm quark were heavy enough that the "inclusive" approach were 
applicable. The hypothesis underlying the present analysis is that this is not the case. 

3 Conclusions 

We proved that if SU(3) violation may be treated perturbatively, then D° — D° mixing in 
the Standard Model is generated only at second order in SU(3) breaking effects. Within 
the exclusive approach, we identified an SU(3) breaking effect, SU(3) violation in final 



6 



state phase space, which can be calculated with minimal model dependence. We found that 
phase space effects alone provide enough SU(3) violation to induce y ~ 10~ 2 . Large effects 
in y appear for decays close to D threshold, where an analytic expansion in SU (3) violation 
is no longer possible. 

Indeed, some degree of cancellation is possible between different multiplets, as would be 
expected in the m c — > oo limit, or between SU{3) breaking in phase space and in matrix 
elements. It is not known how effective these cancellations are, and the most reasonable 
assumption in light of our analysis is that they are not significant enough to result in an 
order of magnitude suppression of y, as they are not enforced by any symmetry arguments. 
Therefore, any future discovery of a D meson width difference should not by itself be 
interpreted as an indication of the breakdown of the Standard Model. 

It is my pleasure to thank S. Bergmann, E. Golowich, Y. Grossman, A. Falk, Z. Ligeti, 
and Y. Nir for collaboration on this and related projects. 



References 

[1] A. Datta, D. Kumbhakar, Z. Phys. C27, 515 (1985); A. Petrov, Phys. Rev. D56, 1685 
(1997). 

[2] A. F. Falk, Y. Nir and A. A. Petrov, JHEP 9912, 019 (1999). 

[3] M. Gronau, Y. Grossman and J. L. Rosner, Phys. Lett. B 508, 37 (2001); J. P. Silva 
and A. Soffer, Phys. Rev. D 61, 112001 (2000). E. Golowich and S. Pakvasa, Phys. 
Lett. B 505, 94 (2001). 

[4] S. Bergmann, Y. Grossman, Z. Ligeti, Y. Nir, A. Petrov, Phys. Lett. B 486, 418 (2000). 

[5] D. Atwood and A. A. Petrov, |arXiv:hep-ph /0207165j 

[6] A. F. Falk, Y. Grossman, Z. Ligeti and A. A. Petrov, Phys. Rev. D 65, 054034 (2002). 

[7] H. N. Nelson, in Proc. of the 19th Intl. Symp. on Photon and Lepton Interactions at 



High Energy LP99 ed. J.A. Jaros and M.E. Peskin, |arXiv:hep-ex/9908021 



[8] H. Georgi, Phys. Lett. B297, 353 (1992); T. Ohl, G. Ricciardi and E. Simmons, Nucl. 
Phys. B403, 605 (1993); I. Bigi and N. Uraltsev, Nucl. Phys. B 592, 92 (2001), for 
a recent review see A. A. Petrov, Proc. of 4th Workshop on Continuous Advances in 
QCD, Minneapolis, Minnesota, 12-14 May 2000, |arXiv:hep-ph/0009160 . 



[9] J. Donoghue, E. Golowich, B. Holstein and J. Trampetic, Phys. Rev. D33, 179 (1986); 
L. Wolfenstein, Phys. Lett. B164, 170 (1985); P. Colangelo, G. Nardulli and N. Paver, 
Phys. Lett. B242, 71 (1990); T.A. Kaeding, Phys. Lett. B357, 151 (1995). 

[10] E. Golowich and A. A. Petrov, Phys. Lett. B 427, 172 (1998). 

[11] D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. C 15, 1 (2000). 



7 



